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Mass loss is a very important aspect of the life of massive stars. After briefly reviewing its 
importance, we discuss the impact of the recently proposed downward revision of mass loss 
rates due to clumping (difficulty to form Wolf-Rayct stars and production of critically rotating 
stars). Although a small reduction might be allowed, large reduction factors around ten are 
disfavoured. 

We then discuss the possibility of significant mass loss at very low metallicity due to stars 
reaching break-up velocities and especially due to the metal enrichment of the surface of the 
star via rotational and convective mixing. This significant mass loss may help the first very 
massive stars avoid the fate of pair-creation supernova, the chemical signature of which is not 
observed in extremely metal poor stars. The chemical composition of the very low metallicity 
winds is very similar to that of the most metal poor star known to date, HE1327-2326 and offer 
an interesting explanation for the origin of the metals in this star. 

We also discuss the importance of mass loss in the context of long and soft gamma-ray bursts 
and pair-creation supernovae. Finally, we would like to stress that mass loss in cooler parts of 
the HR-diagram (luminous blue variable and yellow and red supergiant stages) are much more 
uncertain than in the hot part. More work needs to be done in these areas to better constrain 
the evolution of the most massive stars. 



1 Introduction 

Mass loss has a crucial impact on the evolution of 
massive stars. It affects evolutionary tracks, life- 
times and surface abundances. It also determines 
the population of massive stars (number of stars in 
each Wolf-Rayet subtype for example). It influences 
the type of supernova at the death of the star (SNII, 
lb, Ic, or a pair-creation supernova) and the final 
remnant (neutron star or black hole). Mass loss re- 
leases matter and energy back into the interstellar 
medium in amounts comparable to supernovae (for 
stars above 30 Mq). Finally, it affects the hardness 
of the ionizing radiation coming from massive stars. 
It is therefore very important to understand mass 
loss in order to understand and model the evolution 
of massive stars. 



2 Impact of reduced mass loss 
rates at solar metallicity 

The concept of clumping is not new (see contribution 
from Moffat in this volume, or a general review from 
Kudritzki & Puis 2000). However, new observations 
suggest clumping factors leading to downward revi- 
sion of mass loss rates between three to ten or even 
more for massive stars (Bouret et al. 2005, FuUerton 
et al. 2006). Here we discuss the implications of a 
reduction factor around ten. A 120 Mq star, using 



current mass loss prescriptions at solar metallicity 
(Vink et al. 2000, 2001, Kudritzki & Puis 2000), loses 
on average 2 • 10"^ M© per year. The lifetime of a 
120 Mq star is about 2.5 million years. This implies 
that, on the main sequence, a 120 Mq star loses ap- 
proximately 50 Mq. If mass loss rates are reduced 
by a factor 10, such a star would only lose 5 Mq. 
The question is then, how to produce a WR star with 
such low mass loss rates? The first possibility is that 
mass loss is high in other evolutionary stages, like 
the luminous blue variable (continuum-driven winds 
in LBVs, Smith & Owocki 2006) and the red super- 
giant (RSG) stages. Mass loss rates are harder to 
determine in these two stages and therefore uncer- 
tainties are still large. Nevertheless clumping may 
also affect mass loss determination in other stages. 
Another possibility would be that all massive stars 
are in close binary systems (Kobulnicky et al. 2006). 
However, if this were true, then it would be hard 
to produce the many RSG stars observed. Further- 
more, the fraction of Wolf-Rayet stars in close binary 
systems in the Magellanic Clouds is found to be only 
30-40% (Foellmi et al. 2003, 2003). This means that 
single stars must still be able to lose enough mass 
to become WRs on their own. The last possibility 
discussed here is that rotation (possibly coupled to 
magnetic fields) induces such a strong mixing that 
WRs are produced by mixing rather than mass loss 
(Maeder 1987, Yoon & Langer 2005). This scenario 
works only for fast rotators, which represent only a 
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small fraction at Zq and therefore this cannot pro- 
duce all the WR stars observed. 

Another important impact of strongly reduced 
mass loss rates is that it implies that the angular 
momentum loss is much weaker. This would lead 
to many critically rotating stars near the end of the 
main sequence, similar to Be stars. This in turn 
would lead to an increase in mass loss rates, which 
could possibly compensate for a modest reduction 
factor. Additional models are necessary to give a 
quantitative answer. More interestingly, mass loss 
would become strongly anisotropic (Maeder & Des- 
jacques 2001) and possibly produce disks when the 
fi-limit is reached as is observed around Be type 
stars (See e.g. the review by Porter & Rivinius 2003). 
The lack of observations of critically rotating very 
massive stars is an argument against mass loss being 
extremely low on the main sequence and high only 
later on during the LBV and RSG stages. Note that 
rotation is able to compensate for a reduction factor 
2 for mass loss since comparable number of WR stars 
are produced with enhanced mass loss rates (Meynet 
et al. 1994) and with normal mass loss rates + rota- 
tion (Meynet & Maeder 2005). 



3 Metallicity dependence 

The metallicity (Z) dependence of mass loss rates 
is usually included using the formula: M{Z) = 
M{Zq){Z / Zq)"- . The exponent a varies between 
0.5-0.6 (Kudritzki & Puis 2000, Kudritzki 2002) and 
0.7-0.86 (Vink and collaborators 2001, 2005) for O- 
type and WR stars (See Mokiem et al. 2007 for a re- 
cent comparison between mass loss prescriptions and 
observed mass loss rates). Until very recently, most 
models use at best the total metal content present 
at the surface of the star to determine the mass loss 
rate. However, the surface chemical composition be- 
comes very different from the solar mixture, due ei- 
ther to mass loss in the WR stage or by internal 
mixing (convection and rotation) after the main se- 
quence. It is therefore important to know the con- 
tribution from each chemical species to opacity and 
mass loss. Recent studies (Vink et al. 2000, 2005) 
show that iron is the dominant element concerning 
radiation line-driven mass loss for 0-type and WR 
stars. In the case of WR stars, there is however a 
plateau at low metallicity due to the contributions 
from light elements like carbon, nitrogen and oxy- 
gen (CNO). In the RSG stage, rates generally used 
are still those of Nieuwenhuijzen & de Jager (1990). 
Observations indicate that there is a very weak de- 
pendence of dust-driven mass loss on metallicity and 
that CNO elements and especially nucleation seed 
components like silicon and titanium are dominant 
(Van Loon 2000, 2006, Ferrarotti & Gail 2006). See 
van Loon et al. (2005) for recent mass loss rate pre- 
scriptions in the RSG stage. In particular, the ratio 



of carbon to oxygen is important to determine which 
kind of molecules and dusts form. If the ratio of car- 
bon to oxygen is larger than one, then carbon-rich 
dust would form, and more likely drive a wind since 
they arc more opaque than oxygen-rich dust at low 
metallicity (Hofner & Andersen 2007). 

In between the hot and cool parts of the HR- 
diagram, mass loss is not well understood. Obser- 
vations of the LBV stage indicate that several solar 
masses per year may be lost and there is no indi- 
cation of a metallicity dependence. Chromospheric 
activity could also play a role in stars having sur- 
face temperature similar to the Sun. Thermally 
driven winds and pulsations are still other ways to 
lose mass. Even though there are still large uncer- 
tainties in the dependence of the mass loss rates on 
metallicity in the cooler part of the HR-diagram, it 
is very useful to use models and observations at var- 
ious mctallicities. Indeed, clumping appears to be 
metallicity independent and therefore comparisons 
between models and observations should yield the 
same conclusions at different metallicities. Further- 
more, using models at lower metallicity already give 
a very good estimate of what the impact of clump- 
ing may be on the evolution of the star. The mass 
loss rate is a factor 1.6-2.2 (depending on alpha) and 
2.2-4.0 lower at the metallicity of the large and small 
Magellanic Clouds (LMC and SMC) respectively. 
Comparing models calculated at the SMC metallic- 
ity with observations at solar metallicity shows the 
impact of a reduction factor around three. 

Several groups recently computed massive star 
models and compared them to observed populations 
around solar metallicities. Here we present a few 
of them. Meynet & Maeder (2005) compare the ra- 
tio or WR to 0-type stars using a=0.5 for O-type 
star and no metallicity dependence for WR stars. 
They find that rotating models better reproduce the 
WR/ O ratio and also the ratio of type Ib-|-Ic to type 
II supernova as a function of metallicity compared 
to non-rotating models, which underestimate these 
ratios. Reducing the mass loss rates by even a fac- 
tor two would not fit the observations as well as with 
the current mass loss prescriptions. Eldridge & Vink 
(2006) use mass loss rates dependent on metallic- 
ity in the WR stage and find a better agreement 
with observations for the WC/WN ratio compared 
to metallicity independent mass loss rates in the WR 
stage. Again, reducing the mass loss rates by a fac- 
tor 2 or more would not fit the data better than 
with the current mass loss prescriptions. However 
Vanbeveren et al. (2007) includes binary stars in the 
comparison and find a good fit with a mass loss rate 
reduced by a factor two. 

Including all the arguments discussed above, from 
the current stellar evolution point of view the ob- 
servations of the populations of massive stars would 
not be better reproduced with mass loss rate pre- 
scriptions reduced by a factor greater than two. 
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4 First stellar generations 

As we saw in the previous sections, mass loss plays a 
crucial in the evolution of solar mctallicity stars. In 
this section, we discuss the importance of mass loss 
on the evolution of the first stellar generations. The 
first massive stars died a long time ago and will prob- 
ably never be detected directly (see however Scanna- 
pieco et al. 2005). There are nevertheless indirect ob- 
servational constraints on the first stars coming from 
observations of the most metal poor halo stars (Beers 
& Christlieb 2005). The first stars are very impor- 
tant because they took part in the re-ionisation of 
the universe at the end of the dark ages (roughly 400 
million years after the Big Bang). They are therefore 
tightly linked to the formation of the first structures 
in the universe and can provide valuable information 
about the early evolution of the universe. The first 
stellar generations are different from solar metallic- 
ity (Zq) stars due to their low metal content or ab- 
sence of it. First, very low-Z stars are more com- 
pact due to lower opacity. Second, metal free stars 
burn hydrogen in a core, which is denser and hot- 
ter. This implies that the transition between core 
hydrogen and helium burning is much shorter and 
smoother. Furthermore, hydrogen burns via the pp- 
chain in shell burning. These differences make the 
metal free (first) stars different from the second or 
later generation stars! (Ekstrom & Meynet 2007). 
Third, mass loss is metallicity dependent (at least 
for radiation-driven winds) and therefore mass loss 
is expected to become very small at very low metal- 
licity. Finally, the initial mass function of the first 
steUar generations is expected to be top heavy below 
a certain threshold (Bromm & Loeb 2003). 

Mass loss is expected to be very small. What 
could change this expectation? An additional mech- 
anism or the chemical enrichment of the envelope of 
the star are two possible ways to increase mass loss 
at very low Z. Models of metal free stars including 
the effect of rotation (Ekstrom et al. 2005) show that 
stars may lose up to 10 % of their initial mass due 
to the star rotating at its critical limit (also called 
break-up limit). The mass loss due to the star reach- 
ing the critical limit is non-negligible but at the same 
time not important enough to change drastically the 
fate of the first generation stars. 

The situation is very different at very low but 
non-zero metallicity (Meynet et al. 2006 and Hirschi 
2007). The total mass of an 85 Mq model at Z ^ 
10~^ is shown in Fig. 1 with the top solid line. This 
model, like metal free models, loses around 5% of its 
initial mass when its surface reaches break-up veloc- 
ities in the second part of the main sequence. At the 
end of core H-burning, the core contracts and the en- 
velope expands, thus decreasing the surface velocity 
and its ratio to the critical velocity. The mass loss 
rate becomes very low again until the star crosses 
the HR diagram and reaches the RSG stage. At 
this point the convective envelope dredges up CNO 



elements to the surface increasing its overall mctal- 
licity. The total metallicity, Z, is used in this model 
(including CNO elements) for the metallicity depen- 
dence of the mass loss. 



M= 85 M 
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log|o(Time until core collapse) [yr] 

Figure 1: Structure evolution diagram of the 85 Mq 
model ai Z ~ 10~^. Coloured areas cor- 
responds to convective zones along the la- 
grangian mass coordinate as a function of 
time left until the core collapse. The top 
solid line shows the total mass of the star. 
The burning stage abbreviations are given 
below the time axis. 



Therefore depending on how much CNO is brought 
up to the surface, the mass loss becomes very large 
again. The CNO brought to the surface comes from 
primary C and O produced in He-burning. Rota- 
tional and convective mixing brings these elements 
into the H-burning shell. A large fraction of the C 
and O is then transformed into primary nitrogen via 
the CNO cycle. Additional convective and rotational 
mixing is necessary to bring the primary CNO to the 
surface of the star. The whole process is complex 
and depends on mixing. Multi-dimensional models 
would be very helpful to constrain mixing between 
the hydrogen and carbon rich layers, which releases 
a large amount of energy and strongly affects the 
structure of the star. 

The strongest mass loss occurs in these mod- 
els in the cooler part of the HR diagram. Dust- 
driven winds appear to be metallicity independent 
as long as C-rich dust can form. For this to oc- 
cur, the surface effective temperature needs to be 
low enough (log(Tcff ) < 3.6) and carbon needs to be 
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more abundant than oxygen. Note that nucleation 
seeds (probably involving titanium) are still neces- 
sary to form C-rich dust. It is not clear whether 
extremely low-Z stars will reach such low effective 
temperatures. This depends on the opacity and the 
opacity tables used in these calculations did not ac- 
count for the non-standard mixture of metals (high 
CNO and low iron abundance). It is interesting to 
note that the wind of the 85 Mq model is richer 
in carbon than oxygen, thus allowing C-rich dust 
to form if nucleation seeds are present. There may 
also be other important types of wind, like chromo- 
sphcric activity-driven, pulsation-driven, thermally- 
driven or continuum-driven winds. 

Wind composition, [Fe/H] — 6.6 (Z=10-8) 
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Figure 2: Composition in [X/Fc] of the stellar wind 
for the Z = 10^® models (solid lines). For 
HE1327-2326 {red stars), the best fit for 
the CNO elements is obtained by diluting 
the composition of the wind of the 40 Mq 
model by a factor 600 (see Hirschi 2007 for 
more details). 



Significant mass loss in very low-Z massive stars 
offers an interesting explanation for the strong en- 
richment in CNO elements of the most metal poor 
stars observed in the halo of the galaxy (see Meynet 
et al. 2006 and Hirschi 2007). The most metal poor 
stars known to date, HE1327-2326 (Frebel et al. 
2006) is characterised by very high N, C and O abun- 
dances, high Na, Mg and Al abundances, a weak s- 
process enrichment and depleted lithium. The star 
is not evolved so has not had time to bring self- 
produced CNO elements to its surface and is most 
likely a subgiant. By using one or a few SNe and 



using a very large mass cut, Limongi et al. (2003) 
and Iwamoto et al. (2005) are able to reproduce the 
abundance of most elements. However they are not 
able to reproduce the nitrogen surface abundance of 
HE1327-2326 without rotational mixing. A lot of 
the features of this star are similar to the properties 
of the stellar winds of very metal poor rotating stars. 
HE 1327-2326 could therefore have formed from gas, 
which was mainly enriched by stellar winds of ro- 
tating very low metallicity stars. In this scenario, a 
first generation of stars (PopIII) pollutes the inter- 
stellar medium to very low mctallicities ([Fc/H]~-6). 
Then a PopII.5 star (Hirschi 2005) like the 40 Mq 
model calculated here pollutes (mainly through its 
wind) the interstellar medium out of which HE 1327- 
2326 forms. This would mean that HE1327-2326 is 
a third generation star. In this scenario, the CNO 
abundances are well reproduced, in particular that 
of nitrogen, which according to the new values for a 
subgiant from Frebel et al. (2006) is 0.9 dex higher in 
[X/Fe] than oxygen. This is shown in Fig. 2 where 
the abundances of HE1327-2326 are represented by 
the red stars and the best fit is obtained by diluting 
the composition of the wind of the 40 Mq model by 
a factor 600. When the SN contribution is added, 
the [X/Fe] ratio is usually lower for nitrogen than 
for oxygen. Although the existence of a lower limit 
for the minimum metallicity Z for low mass stars to 
form is still under debate, it is interesting to note 
that the very high CNO yields of the 40 Mq stars 
brings the total metallicity Z above the limit for 
low mass star formation obtained in Bromm & Loeb 
(2003). 



5 Gamma-ray bursts and 
pair-creation supernovae 

Long and soft gamma-ray bursts (GRBs) have now 
been firmly connected to the death of type Ic su- 
pernovae (see Woosley & Bloom 2006 for a recent 
review). In one of the most promising models, the 
collapsar model (Woosley 1993), GRB progenitors 
must form a black hole, lose their hydrogen rich en- 
velope (become a WR) and retain enough angular 
momentum in their core during the pre-supernova 
stages. The strong mass loss discussed in the previ- 
ous section make it possible for single massive stars 
in the first stellar generations to become WR stars 
and even to retain enough angular momentum to 
produce a GRB (Hirschi 2007). A wider grid of 
models at mctallicities around solar shows that the 
rate of fast rotating WO stars is compatible with the 
rate of GRBs with an upper limit around the LMC 
metallicity, in agreement with observations (Hirschi 
et al. 2005). More recent models including the effects 
of magnetic fields (Yoon & Langer 2005) show that 
another mechanism is possible to produce GRBs at 
low Z. This mechanism is the quasi-chemical evo- 
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lution of very fast rotating massive stars. In this 
scenario, WR stars are produced by mixing and not 
mass loss. This last scenario however does not pre- 
dict GRB at metallicities equal or higher than the 
SMC. This upper limit is too low compared to recent 
observations (Fruchter et al. 2006). Taking into ac- 
count the anisotropy in the wind of these fast rotat- 
ing stars (Meynet & Maeder 2007) may help reduce 
the discrepancy between models and observations. 
Note that the downward revision of solar metallic- 
ity (Asplund et al. 2005) may also help resolve the 
problem. 

Apart from GRBs, pair-creation supernovae (PC- 
SNe) are very energetic explosions, which could be 
observed up to very high redshifts (Scannapieco et 
al. 2005). PCSNe are expected to follow the death 
of stars in the mass range between 100 and 250 A'/q, 
assuming that they do not lose a significant frac- 
tion of their mass during the pre-supernova stages. 
Amongst the very first stars formed in the Universe, 
one expects to have PCSN due to the lack of mass 
loss and to the low opacity unable to stop the ac- 
cretion on the star during its formation. However, 
the EMP stars observed in the halo of the galaxy do 
not show the peculiar chemical signature of PCSN 
(strong odd-even effect, see Heger & Woosley 2002). 
This means that either too few or even no PCSN ex- 
isted. One possible explanation to avoid the produc- 
tion of very \ow-Z or metal free PCSNe is the strong 
mass loss in the cool part of the HR diagram due to 
the surface enrichment in CNO elements induced by 
rotational and convective mixing (see previous sec- 
tion) or the star reaching the riF-limit (Ekstrom & 
Meynet 2007). 
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